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Abstract. We report experimental results for heat-transport measurements, in the 
form of the Nusselt number Nu, by turbulent Rayleigh-Benard convection in a 
cylindrical sample of aspect ratio P = D/L = 0.50 {D = 1.12 m is the diameter 
and L = 2.24 m the height). The measurements were made using sulfur hexafluoride 
at pressures up to 19 bars as the fluid. They are for the Rayleigh-number range 
3 X 10^2 < Ra < lO^^^ and for Prandtl numbers Pr between 0.79 and 0.86. For 
Ra < Ral ~ 1.4 x lO^^ gj^^j ^ iVoRa'^-'^^ with -y^ff = 0.312 ± 0.002, consistent 
with classical turbulent Rayleigh-Benard convection in a system with laminar boundary 
layers below the top and above the bottom plate. For Ra^ < Ra < Ra2 (with 
Ra2 ~ 5 X 10^^) 7e// gradually increases up to 0.37 ± 0.01. We argue that above 
Ra2 the system is in the ultimate state of convection where the boundary layers, both 
thermal and kinetic, are also turbulent. 

Several previous measurements for P = 0.50 are re-examined and compared with 
the present results. Some of them show a transition to a state with 7e// in the range 
from 0.37 to 0.40, albeit at values of Ra in the range from 9 x 10^" to 7 x 10^^ which 
is much lower than the present Raj or Ra2. The nature of the transition found by 
them is relatively sharp and does not reveal the wide transition range observed in the 
present work. 

In addition to the results for the genuine Rayleigh-Benard system, we present 
measurements for a sample which was not completely sealed; the small openings 
permitted external currents, imposed by density differences and gravity, to pass 
through the sample. That system showed a sudden decrease of 7e// from 0.308 for 
Ra < Rat ~ 4 X lO^^ ^ o.25 for larger Ra. 

A number of possible experimental effects is examined in a sequence of Appendices; 
none of these effects are found to have a significant influence on the measurements. 
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1. Introduction 

In this paper we consider turbulent convection in a fluid contained between horizontal 
parallel plates and heated from below (Rayleigh-Benard convection or RBC; for a reviews 
written for broad audiences see Refs. [HE]; for more specialized reviews see Refs. |3l|4]). 
The primary purpose of the work on which we report was to search for the transition 
to the "ultimate" state of turbulent convection flrst predicted by Robert Kraichnan [5] 
and Ed Spiegel [6j half a century ago. 

We focus on the particular case of a cylindrical sample of aspect ratio T = D / L = 
0.50 {D = 1.12 m and L = 2.24 m are the diameter and height respectively) because 
this particular geometry was used in previous searches for this state [3 [HI HI [101 [HI [I2l 
[T3t [HI [151 [m [U] ^-iid thus enables a more direct comparison with earlier measurements. 
Experiments searching for the ultimate state using other values of F are of course 
important as well and some have been carried out |T71 UHl [191 [201 [211 [22]; but they 
are beyond the scope of this paper. The work reported here consists of measurements of 
the heat transport by the turbulent system. Other aspects will be discussed separately. 

We present results that were obtained in the High-Pressure Convection Facility 
(the HPCF, a cylindrical sample of 1.12 m diameter) at the Max Planck Institute for 
Dynamics and Self-organization in Gottingen, Germany using sulfur hexafluoride (SFg) 
at pressures up to 19 bars as the fluid. Early results from this work were presented in 
Refs. [231 [Ml [25]. A description of the apparatus was given in Ref. [23]. The present 
paper presents new results obtained after various sample-chamber modiflcations to be 
described in this paper and is a comprehensive report on this work. A brief report of 
these recent results was provided in Ref. [26] . 

The HPCF is located inside a pressure vessel known as the Uboot of Gottingen 
which gets fllled with the gas of choice. Thus there must be some way for the gas to 
enter or leave the HPCF. Originally a small gap of average width about one mm was 
permitted for this purpose to remain between the top and bottom plate and the side 
wall of the sample [23] . Erroneously it was assumed that this gap was negligible because 
it is small compared to the 1120 mm sample diameter. This sample will be called the 
"open" sample. It turned out that the results depended on the temperature difference 
between the fluid in the sample at a temperature Tm and the fluid in the remainder of 
the Uboot at Tf/ [27]. For that reason the HPCF was modifled by sealing the bottom 
plate to the side wall, but leaving the gap between the top plate and the side wall open. 
This system will be called the "half-open" sample. There still was a major difference 
between the results for Tm > Tu and Tm < Tu, although these results differed from 
those of the closed sample. Thus, as a flnal measure, both the top and the bottom 
plates were sealed to the side wall, and a 25 mm diameter tube was installed to permit 
the gas to enter the HPCF. One end of the tube was flush with the inside of the side 
wall, and the other terminated in a remotely controlled valve. The sample could thus 
be fllled while the valve was open, and then during measurements the valve could be 
closed. This sample will be called the "closed" sample. Results for all three versions 
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will be presented in this paper. 

In the next section we shall define the parameters needed to describe this system 
(Sec. E]). We shall then, in Sec. [31 outhne the main features of turbulent convection as 
they are now understood. First, in Sec. 13. we describe the classical state of turbulent 
convection which exists below the transition to the ultimate state with turbulent BLs. 
This will be followed in Sec. 13.21 by a description of what is known or expected for the 
ultimate state. This introductory material will be followed in Sec. |4]by a brief discussion 
of the apparatus modifications used in this work. A detailed description of the main 
features was presented before [23]. Section [5] presents a comprehensive discussion of our 
results and of the results of others at large Ra for cylindrical samples with F = 0.50. It 
is followed by a Summary of our work in Sec. [61 

In a sequence of six Appendices we discuss a number of experimental issues which 
might be of lesser interest to the general reader but which are of considerable importance 



to the specialist. First, in Appendix A we compare measurements for three different 
angles of tilt of the sample axis relative to gravity. The effect of a tilt was studied before 
by several groups, with varying results [28l [29l |30l [311 [321 117]. A tilt is used at times by 
experimentalists to give the large-scale circulation in the sample a preferred azimuthal 
orientation. This was our motivation as well. We show that it had no discernible effect 



on Nu. In Appendix B we demonstrate that non-Oberbeck-Boussinesq effects [331 E 



have only a very minor effect on Nu which can be seen only at the largest values of AT 



near 20 K. In Appendix C we present values of the parameter ^ introduced recently by 
Niemela and Sreenivasan ^35j to describe a special non-Oberbeck-Boussinesq effect which 
apparently occurs near critical points. We find that for our work ^ ~ 1.3, indicating that 
the effect discussed in Ref. [35] (which occurs for small or negative ^) is not expected 



to be relevant to the present work. In Appendix D we give results for the horizontal 



temperature variations in the top and bottom plate of the sample and demonstrate 



that they do not influence the measured values of Nu. In Appendix E we show that 



the "closed" sample really was completely sealed. In Appendix F we present data for 
the influence of a mismatch between the temperature of the sample side shield and the 
sample mean temperature T^, and show that the side-shield temperature-regulation 



is good enough for these effects to be negligible. Finally, in [Appendix G[ we give a 
complete list of our data in numerical form. 

2. The system parameters 

For turbulent RBC in cylindrical containers there are two parameters which, in addition 
to F, are expected to determine its state. They are the dimensionless temperature 
difference as expressed by the Rayleigh number 

Ra.^^ (1) 

and the ratio of viscous to thermal dissipation as given by the Prandtl number 

Pr = u/k . (2) 
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Here a is the isobaric thermal expansion coefficient, g the gravitational acceleration, 
K the thermal diffusivity, u the kinematic viscosity, and AT = — Tf the applied 
temperature difference between the bottom (T^) and the top (Tt) plate. 

In the present paper we present measurements of the heat transport. These results 
are presented in the form of the scaled effective thermal conductivity known as the 
Nusselt number which is given by 

Nu ^ . (3) 

Here Q is the applied heat current, A = D'^tt/A the sample cross-sectional area, and A 
the thermal conductivity. The measurements cover the range 10^^ < Ra < 10^^ and are 
for Pr ranging from 0.79 at the lowest to 0.86 at the highest Ra. 

All fluid properties needed to calculate Ra, Pr, and Nu were evaluated at the mean 
temperature = (Tt + Tf,)/2 of the sample. They were obtained from numerous papers 
in the hterature, as discussed in Ref. I36l. 



3. The characteristics of turbulent RBC 



3.1. The classical state. 

A "classical" state of RBC exists below a transition range to an "ultimate" state; the 
transition range extends over more than a decade, approximately from Raj to Ra2 p6] . 
For simplicity of discussion we shall characterize this range by Ra*, taken to lie perhaps 
somewhere near the middle of the range, which, for the parameters of our work, is about 
10^^ inZl For Ra < Ra* the heat transport in this system is controlled by laminar 
thermal boundary layers (BLs), one just below the top and the other just above the 
bottom plate. The value of Ra* has been the subject of discussion for some time, and a 
major issue at the forefront of the field is the nature of the state above Ra*. Estimates 
of Ra* are not very accurate; a reasonable argument |37] yielded Ra* ~ 10^^ or so for 
Pr ~ 1, although another estimate [19] gave a value closer to 10^^. 

For Ra < Ra* nearly half of AT is found across each BL, and the sample interior 
(known as the "bulk") has a highly fluctuating temperature which is nearly uniform 
in the time average [381 EH]. At a more detailed level it was recognized long ago that 
the bulk actually sustains small temperature gradients, but the total temperature drop 
across it is much smaller than that across the BLs (see, for instance, [iQl UH Il2]). Very 
recently it was found that these small temperature variations in the bulk take the form 
of a logarithmic dependence on the distance from the plates [13]; but the precise origin 
of this logarithmic variation is not yet known. 

For the classical state it is well established both experimentally [HJ SHI HHl [El HTl 
SH 119] and theoretically [50] [511 |37l [52] that the Nusselt number can be represented by 
a power law 

Nu = iVoRa^-^^^ (4) 
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with the effective exponent varying from about 0.28 near Ra = 10^ to about 0.32 near 
Ra = 10^^, at least when Pr is close to one or larger. 

It is also well established that, in cylinders with F ~ 0.5 containing a fluid with 
Pr ~ 0.7 and for Ra < 10^^, there is a large-scale circulation (LSC) in the sample 
interior that takes the form of a single convection roll, with up-flow along the side wall 
at an azimuthal position 6o and down-flow also along the wall but at an azimuthal 
position close to + ^ ISHl 112] • The LSC is bombarded by the small-scale fluctuations 
of the system, and may be regarded as a stochastically driven system that fluctuates 
intensely [5l| 155] . Even at modest Ra below, say, 10^^ and for Pr near 0.7, the LSC 
existence is intermittent. It frequently collapses, only to re-form again at a somewhat 
later time. Whether the LSC survives at all up to Ra* had not been clear heretofore; 
it was recognized (see, for instance, [19j) that the LSC becomes less well defined as Ra 
increases, but concrete quantitative experimental evidence for its existence or demise 
has only become available during our present work. We found that, even for Ra = 10^^, 
there is evidence of its existence, but its average lifetime is short and it may be regarded 
more appropriately as just one of the continuum of modes contributing to the fluctuation 
spectrum of the system. This will be reported in detail in a subsequent paper. 

When the fluctuations are not too vigorous, the LSC due to its horizontal flow with 
speed U just adjacent to the top and bottom plate will establish viscous BLs adjacent 
to the plates. The viscous BLs may be imbedded in the thermal ones or vice versa, 
depending on Pr. These BLs are laminar, albeit fluctuating [561 EZ], in the classical 
state. 

3.2. The ultimate state. 

About half a century ago it was predicted by Kraichnan [5| and Spiegel ^ that, in the 
absence of boundary layers, the Nusselt number should be proportional to Ra'^ with 
7 = 1/2. This prediction is consistent with rigorous upper bounds for Nu(Ra) obtained 
by Howard [58j and by Doering and Constantin [5^. Although it seems difficult to 
construct a physical system without boundaries, the 7 = 1/2 prediction was supported 
by direct numerical simulations (DNS) of RBC with periodic boundary conditions (BCs) 
in the vertical direction and forcing in the bulk [601 [611 162] . as well as by DNS for the 
Rayleigh- Taylor instability [63] which is expected to reveal similar phenomena and has 
no boundaries. Experimentally, it is noteworthy that a local heat-flux measurement in 
the center of a Rayleigh- Benard sample in the classical state {i.e. in the state with 
laminar BLs) yielded an exponent of 0.5 [63], even though the global heat flux led to 
7e// ^ 0.3. 

In the presence of boundaries, Kraichnan noted that the BLs should become 
turbulent when Ra exceeds some characteristic value Ra*. This event was expected 
to be induced by the shear applied to the BLs by the large-scale circulation, or if none 
exists, by the vigorous turbulent fluctuations in the sample interior. It was expected to 
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occur when the shear Reynolds number 

Re, = UXJu (5) 

determined by U and the BL thicknesses A„ exceeded a typical value somewhere in 
the range from about 200 to 400 [65]. When there is no LSC because fluctuations 
dominate, one would expect the fluctuations to take on its role and generate shear, 
with their root-mean-square velocity V near the BLs taking the role of U in Eq. |5l 
In that case the characteristic size of fluctuations will cover a range, roughly from D 
down to smaller lengths, and will be intermittent in time. One then would expect 
the turbulent shear layers to be more localized laterally in space, as well as in time. 
In Kraichnan's considerations he assumed that the viscous and thermal BLs would 
undergo the shear-induced transition at the same value of Ra. Even in the presence of 
rigid top and bottom plates the prediction for the large-Ra asymptotic state then was 
still Nu ~ Ra^^^. However, Kraichnan [5] predicted that, due to the turbulent BLs, 
there would be logarithmic corrections to this power law. 

Recently Grossmann and Lohse considered the possibility that the thermal and 
viscous BLs may undergo the turbulent shear transition at different values of Ra or 
simultaneously, and derived the consequences of transitions in one or the other or both 
[66] . The ultimate state would then correspond to the case where both the viscous and 
the thermal BLs become turbulent (see Sec. III.C of [66]; we will reserve the notation 
Ra* for this case). 

In analogy to the logarithmic velocity profiles in turbulent shear flows first 
considered by von Karman [67] and Prandtl [68] (for a recent review, see [69]). 
Grossmann and Lohse [66] predicted that the turbulent BLs would extend throughout 
the sample, replacing the bulk by a temperature profile that varies logarithmically with 
the distance from the plates. For the Boussinesq system [331 El] the two profiles, one 
coming from the top and the other from the bottom plate, would then meet at the 
horizontal mid-plane of the cell. Logarithmic temperature profiles have indeed been 
observed in recent measurements for the ultimate state [13]; but since they were found 
for the classical state as well, it remains unclear to what extend this finding supports the 
prediction. In conjunction with the viscous and thermal sublayers near the plates (which 
survive above Ra* because of the boundary conditions at the solid-fluid interface) the 
extended turbulent BLs lead to logarithmic corrections to the asymptotic power law for 
Nu. The Grossmann and Lohse prediction for these logarithmic corrections differs from 
the original form of the logarithms given by Kraichnan. However, for either prediction 
the corrections vary only slowly with Ra, and in experimentally accessible Ra ranges 
one expects an effective power law with an effective exponent 7e// — 0.38 to 0.40. The 
asymptotic regime where the effective exponent has essentially reached 1/2 is well out 
of reach of any conceivable experiment. 

An important question was the value of Ra*. Since the Reynolds number Re = 
UL/v of the global LSC (and thus U and Re, in Eq. |5]) decrease with increasing Pr at 
constant Ra (see, for instance, Ref. [37]) it follows that Ra* increases with increasing 
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Pr. The value of Re* depends on the nature and amphtude of prevaihng perturbations, 
but is estimated to be in the range from 200 to 400 [65] . For Pr ^ 1 these considerations 
led to Ra* = C(10^^) [37]; but other estimates [I9] yielded lower values. 

A notable recent success in the search for the ultimate state has been achieved 
with turbulent Couette- Taylor (CT) flow [TOl [71] (it had been shown on the basis of its 
equations of motion that it should undergo an ultimate-state transition that is analogous 
to that of RBC [72]). In the CT case the shear is applied directly to the fluid by 
concentric rotating cylinders, and thus is much more effective in driving the BLs into 
the turbulent state than is the shear in RBC which is generated as a secondary effect 
by the buoyancy-induced LSC or the fluctuations. The CT measurements yielded an 
effective exponent of 0.38 for the corresponding variables, consistent with 1/2 and the 
predicted logarithmic corrections [66] . 

For RBC the situation is less clear. In order to reach exceptionally high Ra, 
two groups used fluid helium near its critical point at temperatures of about 5 K and 
pressures of about 2 bars. One of them [HI [12], at the time located at Grenoble, reported 
to have found the ultimate regime, and cited a value Ra* ~ 10^^. We shall refer to these 
results as the "Grenoble" data. A major puzzle created by these results is that one can 
estimate that the data imply Re* ~ 100 or less, and this seems too low for any BL shear 
instability. 

In a second nominally equivalent investigation near the critical point of helium 
Niemela et al. [10] made measurements of Nu up to Ra ~ 10^^. They found that 
Nu oc Ra*^'^^ [73] up to their largest Ra, without any evidence for a transition. This 
work was done at the University of Oregon, and we shall refer to these results as the 
"Oregon" data. In this case the absence of a transition does not necessarily contradict 
expectations because Pr began to increase as Ra exceeded about 10^^, and it is plausible 
that Ra*(Pr) never was reached or resolved in that experiment. 

There has been a number of additional low-temperature experiments intended to 
clarify the situation; we refer to a recent review [3] for a detailed discussion of these 
measurements. For completeness we mention a comprehensive recent article by Roche 
et al. [17] which examines the influence of the nature of the side walls, of F, of Pr, 
and of several other factors which seem to influence the transition to a state with 
7e// significantly larger than 0.32. This survey concludes that a transition to the 
ultimate regime occurs in several experiments, but again these transitions occurred 
at unexpectedly low values of Ra. The survey concludes that the transition occurs at 
smaller Ra when Pr is larger, which is opposite to the Pr dependence of Ra* expected for 
the shear instability. Since no LSC-induced shear instability is likely to have occurred, 
and since the ultimate-state predictions are based on the assumption of turbulent BLs, it 
remains unclear to us how the states with 7e// much larger than 0.32 reported in Ref. [17] 
are related to the Kraichnan prediction [5] or to the states discussed by Grossmann and 
Lohse [66l. 
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Table 1. Versions of the High-Pressure-Convection Facility (HPCF). The second 
column gives the material used to construct the top and bottom plates. The third 
column indicates whether the top and/or bottom plate was sealed to the side wall. 



Version 



Plates 



Seals 



runs 



HPCF-I 
HPCF-IIa 
HPCF-IIb 
HPCF-IIc 
HPCF-IId 
HPCF-IIe 



Aluminum 



Copper 
Copper 
Copper 
Copper 
Copper 



bottom and top 100918 to 110919 



bottom and top 100612 to 100818 



none 080827 to 090313 

none 090505 to 090917 

none 090905 to 100125 

bottom 100202 to 100502 



4. Apparatus 

Versions HPCF-I, HPCF-IIa, and HPCF-IIb of the apparatus were described in detail in 
Ref. [23]. A schematic diagram of these units can be found in Fig. 2 of that reference. 
HPCF-I had aluminum top and bottom plates, whereas HPCF-II had copper plates. 
HPCF-IIa and lib differed only in the amount and type of insulation (mostly open-pore 
foam and aluminum-coated polyester film, see p3]) provided outside the sample cell, and 
data obtained with them showed no obvious difference. More recently three additional 
modifications known as HPCF-IIc, -Ild, and -He were developed, and corresponding 
schematic diagrams are shown in Fig. [H HPCF-IId and -He differed from HPCF-IIc 
only by the 2.5 cm diameter side arm and valve used to fill and empty the sealed samples, 
and by whether or not the side wall was sealed to the bottom and/or top plate. Table [T] 
lists the major differences between the six versions, as well as the identifications of the 
run^ performed in each. 

All samples had an internal height L = 2240±2 mm and diameter D = 1122±2 mm. 
The measurements to be discussed here were made in HPCF-IIb to -He. All samples had 
a Plexiglas side wall of 0.95 cm thickness and several thermal shields. The entire sample 
was immersed in a high-pressure vessel, known as the Uboot of Gottingen, that could 
be filled with various gases, including sulfur hexa- fluoride (SFg), up to a pressure of 
P = 19 bars. As shown in Fig. [1], all samples had a composite bottom plate consisting 
of a bottom (BPb) and a top (BPt) member made of copper (aluminum for HPCF- 
I) and a 5 mm thick layer of Lexan sandwiched between them. The composite was 
glued together with very thin layers of degassed Stycast 1266 epoxy. The temperature 
difference across this composite, together with the composite conductance, was used 
to infer the heat current Q that entered the sample at the bottom. The underside of 
the bottom member of the composite was heated electrically by a heater immersed in 
epoxy in grooves ("BP htr grooves"). The top plate (TP) was cooled by a water circuit 
consisting of two pairs of spirals. The pairs were in parallel, and the flows in the two 

I The run numbers had the structure "yymmdd" . 
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Figure 1. Left: Diagram of HPCF-IIe. Right: Detailed diagram of the top and bottom 
corners of HPCF-IIc to HPCF-IIe. All parts are shown to scale, except for the valve 
in the left part. We refer to Fig. 2 of 23 for a description of many features that were 
common with HPCF-I to HPCF-IIb. The bottom plate was a composite consisting 
of a bottom member "BPb" , a top member "BPt" , and a Lexan sheet between them. 
The bottom-plate heater was imbedded in epoxy inside the heater grooves ("BP htr 
grooves"). The bottom shield ("BS") was extended by adding a section ("BS2"). The 
bottom bulk microshield ("BBMS"), servoed at Tm, is new. The side shield ("SS") 
and top microshield ("TMS") are unchanged. The location and size of the spiral water 
channels ( "water spirals" ) in the top plate ( "TP" ) are indicated. 

members of a given pair were anti-parallel. Remaining horizontal thermal gradients in 
the TP and the BPt are discussed in Appendix D. 

The various shields which prevented parasitic heat losses from the sample cell were 
discussed in detail in [23]. Starting with HPCF-IIc we added two more shields, the 
bottom shield, which is always servoed at the temperature of the bottom member BPb 
of the composite, was extended by adding the section BS2 (see Fig.[T]). A more significant 
addition was the "Bottom Bulk Micro Shield" BBMS. It was servoed at Tm and thus 
minimized vertical thermal gradients in the space between the side wall and the side 
shield (SS). Prior to the addition of BBMS there was a vertical temperature drop from 
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the BPt temperature at BMSt to at TMS, which is approximately equal to AT/2 
and thus generated a Rayleigh number about equal to half the sample Ra. Even though 
the space between the side wall and the SS was filled with foam and polyester film, 
convection was believed to have been induced in this space during runs at the larger 
values of Ra. 

HPCF-IIc was identical to HPCF-IIb, except that a seal consisting of silicone 
adhesive was applied to the inside corner between the side wall and the top of the 
bottom-plate composite along the entire periphery. It is expected that this seal will 
prevent any fiow through the small gap, of width about 1 mm, between the side wall 
and the bottom plate. A similar gap between the side wall and the top plate was left 
open since fiuid had to be allowed to enter or leave the cell as the temperature or pressure 
was changed. 

HPCF-IId and HPCF-IIe consisted of a completely sealed system, with no gaps 
between the top or bottom plates and the side wall. A tube of inside diameter 2.5 
cm was installed and entered the side wall at half height. Its termination was fiush 
with the inside of the side wall, without any protrusion into the convection chamber. 
Outside the convection chamber this tube contained a remotely operated ball valve. A 
small-diameter (~ 3 mm) tube led from the sample side of the 2.5 cm diameter tube to 
a location outside the Uboot where it was connected to the pressure gage. Thus, the 
actual sample pressure could be monitored. At each set point of the experiment the 
system was equilibrated with the valve open for about six hours. The valve was then 
closed, permitting measurements on a completely sealed system. The fill tube had two 
side arms with an additional valve in each. One opened when the pressure difference 
Pjj — P between the Uboot and the sample exceeded 25 mBar; the other opened when 
this pressure difference was less than -25 mBar. We note that a pressure difference of 
25 mBar leads to a force of about 250 kg acting on the top and bottom plates. When 
for instance this force exceeds the weight of the top plate, then this plate will lift up 
and damage will be done to the instrument. In order to keep \Pu — P\ sufficiently small, 
filling and emptying of the Uboot and sample was done very slowly, over a period of 
a day or two depending on the desired pressure. Measurements of the sample pressure 
under various conditions showed that the sample was indeed sealed, as discussed in 
Appendix E. 

Nusselt-number measurements were based on temperatures determined with fifteen 
thermometers, five each in the two members BPb and BPt of the bottom-plate composite 
and five in the top plate TP. Each set of five consisted of one thermometer placed at 
the plate center and four, positioned equally spaced azimuthally, at a distance of 0.42£) 
from the center. The thermometers in the BPt and the TP were located about 1 mm 
from the fiuid-copper interface. The three sets of five thermometers were used to obtain 
the averages Tspb, Tspt, and T^p. Tppb and Tppt were used to calculate the heat current 
entering the sample. Tspt and T^p were used to obtain AT and Ra. In a typical run 
both Tppb and Ttp were regulated at a specified setpoint. 

A small correction to Nu was made for the sidewall conductance [TU [75]. This 
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Figure 2. The Nusselt number Nu as a function of the Rayleigh number Ra on 
logarithmic scales. Solid black circles: This work, T„i — Tu< — 2K. Open black circles: 
This work, Tm-Tu > 2K. Red stars: the Oregon data [TOl El [73] . Open blue circles: 
the Grenoble data Blue sohd line: the power law Nu = 0.1044Ra°-^^^. 



correction was about 1.4% for Ra ^ 10^^ and about 3.5% for Ra ~ 5 x 10^^. Neglecting 
this correction changed the exponent obtained from a power-law fit to the data for 
Nu(Ra) by about -0.003. Estimates and a comparison of measurements with copper 
and aluminum plates [76] showed that corrections for the finite conductivity of the top 
and bottom plates [771 EH] were negligible. 

5. Results. 

5.1. Closed sample. 

5.1.1. The broad overview. In this subsection we present results for the closed sample 
HPCF-IIe which is our best approximation to the idealized RB system. In Sect. 15.21 
we compare these results with previous measurements. Then, in Sec. 15.31 "we discuss 
the measurements for the open sample where the RB system may be perturbed by an 
additional current entering or leaving the sample through the narrow gap between the 
top and bottom plates and the side wall because of the chimney effect. Finally, in 
Sec. 15.41 we consider the case where only the bottom of the sample is sealed (HPCF-IIc) 
while a gap remains between the top plate and the side wall (the "half-open" sample). 

Results of the Nu measurements for the closed sample (HPCF-IIe) were reported 
briefly before [26]. It was found that they depended slightly on Tm — Tu, but much 
less so than the data for the open or the half-open sample. We have been unable to 
determine the reason for this dependence which persisted in spite of the many thermal 
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Figure 3. a): The reduced Nusselt number Nu/Ra^'^"'^^ on a linear scale as a function 
of Ra on a logarithmic scale. Solid (open) black circles: T,„ — Tij< —2 K (T^ — Tjj > 2 
K). The points indicated by a red dot in the center were used in a least-squares fit of 
a power law to the data. 



shields and the foam and foil insulation that were provided (see [23], and Sec. [Hand 
Fig. [1] above). 

In Fig. [2] we show Nu as a function of Ra with both axes on a logarithmic scale. 
The open circles are for — T^/ > 2 K, and the solid black circles represent the data 
for Tm — Tu < —2 K. Within the resolution of this figure the open and solid circles are 
seen to agree quite well with each other and with the Oregon data (red stars), although 
small differences can be noticed on close inspection. Both our data sets and the Oregon 
data differ significantly from the Grenoble data (open blue circles). A more detailed 
comparison with those as well as with other [79l |17] results is given below. 

5.1.2. The classical state. First we examine Nu(Ra) in the classical regime in greater 
detail. Figure [3] is a high-resolution graph of the data in this parameter range in the 
form of Nu/Ra°'^^^ on a linear scale as a function of Ra on a logarithmic scale. One 
sees that each of the two data sets covers a range of about a decade in the classical 
regime where a simple power law describes them well. A fit of the power law Eq. H] to 
the data points that are indicated by small red dots in their centers gave the parameter 
values and standard erroers (67% confidence limits) listed in Table El On the basis of 
this analysis we chose 7e/j = 0.312 as our best estimate of the effective exponent in the 
classical regime. As can be seen from the table, the uncertainty of this result due to 
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the scatter in the data is less than 0.001. We estimate that an additional uncertainty 
comes from possible systematic errors of the sidewall correction, and thus the overall 
uncertainty of the exponent is about 0.002. This value is consistent with numerous other 
measurements at smaller Ra and larger Pr, and agrees quite well with the value 0.323 
obtained from a numerical analysis of the Grossmann-Lohse model, Eqs. 13 and 14 of 
Ref. [51], for 10^^ < Ra < lO^^ and Pr = 0.8 (but for F = 1.00). 

Table 2. Results of least-squares fits of a power law to the data in the classical regime. 

Data set Nq •y^ff 

Tm-Tu< -2K 0.1040 ±0.0011 0.3121 ± 0.0004 

Tm-Tu< -2K 0.1044 ± 0.00002 0.312 (fixed) 

T^-Tu > 2K 0.1020 ±0.0037 0.3116 ± 0.0012 

Tm-Tu> 2K 0.1006 ± 0.00009 0.312 (fixed) 



5.1.3. Transition to the ultimate state. In order to explore the dependence of the data 
on Tm — Tu in more detail we show in Fig. the results for the reduced Nusselt number 
Nured = Nu/ Ra^'^^"^ on a linear scale as a function of Ra on a logarithmic scale. Here 
it becomes apparent that the Tm — Tu < —2 K data (solid circles) are higher than the 
Tm — Tu > 2 K data (open circles) by about 6% near Ra = 10^"^ and about 10% near 
Ra = 10^^. In the classical regime Ra < 10^'^ (see Fig. [3] and Table [2]) the difference is 
3.8%. In Fig. nil we added a third set of data taken at nearly constant AT ~ 10.3 K but 
at various values of — Tu- During these measurements Tu was not controlled and 
determined by the balance between the heat input to the Uboot from the HPCF-II and 
the cooling to the surrounding laboratory. It varied over the narrow range from 24.3 to 
25.5°C. The sample temperature was controlled by feedback loops and was changed 
in small steps from 21 to 27°C. Since Tm (and thus the fiuid properties) changed, the 
data at constant AT led to a small variation of Ra. One sees that they cover the Nu^ed 
range from the upper to the lower branch. 

The same constant AT results are shown also in Fig. |Dd, but as a function of T^—Tu- 
Here one sees that the data become independent of Tm — Tu when \Tm — Tu\ > 2K. This 
is the reason why the majority of data (the open and solid black circles) were taken as 
a function of Ra with \Tm — Tu\ > 2 K. 

In Fig. [5] we show all the data for the two states with \Tm — Tu\ > 2 K over 
the entire Ra range accessible in the experiment. One sees the classical state with 
Jeff = 0.312 for Ra < 10^'^. At larger Ra the two data sets trace out curves with 
remarkably similar shapes, albeit displaced both vertically and horizontally. To further 
explore the similarity between the two sets we multiplied the open circles by 1.04 . 
This yielded the open blue squares, which now agree with the solid black circles in the 
classical range. Further dividing the Ra values of the open blue squares by 3.7 yielded 
the open red diamonds. One sees that these two transformations yielded agreement 
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Figure 4. a): The reduced Nusselt number Nu/Ra° '^^^ on a linear scale as a function 
of Ra on a logarithmic scale. Solid (open) black circles: T,„ — T[/ < —2 K {Tm — Tjj > 2 
K). Sohd red squares: Nu,.ed = Nu/Ra°-^^^ at nearly constant AT ~ 10.3 K. b): 
Nu/Ra°'^^^ at nearly constant AT ~ 10.3 K as a function of T,„ — T[/ on linear scales. 



within the experimental scatter between the data at large and small Tm — Tjj. This 
shows that the shapes of the curves traced out by the two data sets are the same. 

Both data sets reveal a departure from the classical effective power law, with Nu 
increasing more rapidly with Ra than Ra°'^^^ when Ra > Raj where Raj ^ 1.5 x 10^^ 
for the solid circles and ~ 5 x 10^^ for the open circles. Henceforth we shall concentrate 
on the results for Tm — Tu<—2 K. They continue to increase beyond the classical-state 
values, with an effective exponent that gradually becomes larger until Raj ~ 5 x 10^^ is 
reached. Beyond Rag one has ■jeff — 0.37 ± 0.01 as shown by the blue solid line in the 
figure. This result is consistent with the prediction of an asymptotic exponent 7 = 1/2 
modified by logarithmic corrections in the ultimate state with turbulent boundary layers 
above the bottom and below the top plate. The recent prediction by Grossmann 
and Lohse [SS] for the form of the logarithmic corrections differs from that given by 
Kraichnan [5]; but our data can not distinguish between these two theoretical results 
which both yield values of 7e// which are roughly in the range from 0.38 to 0.40. 

It is worth noting that the data in the transition range Raj < Ra < Rag have 
significantly greater scatter than the data in the classical regime Ra < Raj or those 
in the ultimate regime Ra > Rag. This indicates the existence of multiple states, 
presumably with subtly different BL configurations, during the complex transition from 
laminar to turbulent BLs. 

We call attention to the fact that the transition range between the classical and 
the ultimate state can also be found, between about the same values of Raj and Rag, 
in results for the Reynolds number and in measurements of vertical logarithmic 
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Figure 5. The reduced Nusselt number Nu/Ra" '^^^ as a function of Ra for tlie "closed" 
sample. Black solid circles: Tm — Tij < — 2K. Black open circles: T,n~Tu > +2K. Open 
squares (blue online): Nu values of the open black circles multiplied by 1.04. Open 
diamonds (red online): Ra values of the open squares divided by 3.7. The horizontal 
lines represent power laws for Nu with j^ff = 0.312 and A^o = 0.1005 and 0.1044. The 
solid (dotted) line (blue online) through (near) the data at the largest Ra corresponds 
to 7e// = 0.37. 

temperature profiles that extend over much of the sample height [43]. 
5.2. Comparison with previous results. 

5.2.1. Comparison with the Crenohle, Oregon, and DNS data. Throughout this 
comparison we shall focus on our data for — Tu < — 2K; they will be shown as 
black solid circles in Figs. [6] and El We mention again that, for our results, Pr changed 
only over the narrow range from 0.79 near Ra = 3 x 10^^ to 0.86 near Ra = 10^^. 
Throughout this comparison we shall show the data of others in red for Pr < 1.0, green 
for 1 < Pr < 2, blue for 2 < Pr < 4, and purple for 4 < Pr < 8. 

In Fig. [6] our measurements are compared with those of Niemela et al. [TD] (Oregon 
data) which are given as stars. For the large range 10^ < Ra < 2 x 10^^ those 
results agree well with the power-law fit to our data in the classical regime which gave 
Nu = 0.404Ra°-^^l One can argue that they rise above this power law as Ra approaches 
10^^, perhaps signaling the beginning of a transition to the ultimate state, but in view 
of the scatter of the data this argument might not be convincing. It is noteworthy that 
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Figure 6. Comparison of the present results (solid black circles) with the Oregon 
(stars) [TU], the Grenoble [T2] (small open symbols), and the direct numerical 
simulation fSIT (DNS) data. For the Oregon and Grenoble data we used different 
colors for different ranges of Pr. Red: Pr < 1. Green: 1 < Pr < 2. Blue: 2 < Pr < 4. 
Purple: 4 < Pr < 8. The DNS data are for Pr = 0.7 and are shown as purple circles 
with plusses and error bars. The slanting green and blue line are power-law fits to the 
points with small dots in their centers. 

the rise occurs when Pr is still less that one. However, as Ra increases beyond 10^'^, no 
further increase above the classical power law occurs. It may be that this saturation 
is due to an increase of Pr, which begins to occur just in this Ra range and which is 
expected to shift Ra* to higher values. In summary, the Oregon data show a departure 
from the classical power law near Ra = 10^^, but the evidence for having entered the 
ultimate state in our view remains inconclusive. We note that the original authors of this 
work believed that their data could be represented within their scatter and accuracy by 
a single power law with 7e// = 0.32 over the entire Ra range up to 10^'', thus providing 
no evidence for an ultimate-state transition. 

Also shown in Fig. [6l as small open symbols, are the Grenoble data [12j that fall in 
the range of the figure. One can see that they form distinct groups, depending on Pr. 
For 1 < Pr < 2 (green symbols) there are data in the classical regime. They yield the 
effective exponent of 0.317 ± 0.003, which is consistent with 7e// = 0.312 as adopted 
by us when possible systematic errors, particularly due to uncertainties in the side-wall 
corrections |74l |75], are taken into consideration. The actual values of Nq = Nu/Ra°'^^^ 
are also quite close to our result of 0.104. An average of the data for 10^ < Ra < 10^^ 
is only 2% higher, as shown by the horizontal green line in the figure. 

For 1 < Pr < 2 there is a sharp transition to a Ra range where the effective exponent 
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Figure 7. (a): A higher- resolution graph, on double logarithmic scales, of the 
Grenoble data [IJI for Pr = 1.1 and 1.3. A power-law fit to the blue circles (red 
diamonds) yielded 7e// = 0.317 ± 0.003 (7e// = 0.382 ± 0.006). The green squares 
were not used in either fit. (b): Our data are plotted with the same resolution as that 
used in (a), but with the horizontal axis shifted by two decades. A power-law fit to 
the blue circles (red diamonds) yielded 7e// = 0.3121 ± 0.0004 (7^// = 0.371 ± 0.01). 
Our data show a transition region with a width of about 1.5 decades (green squares), 
whereas the Grenoble data suggest a sharp transition. 



is larger than the classical value of 0.312. In order to determine the Rayleigh number Rat 
at the transition and the effective exponent above it objectively, we selected a subset of 
points which we deemed to be above Ka-t and which seemed consistent with an effective 
power law. These points are indicated by a small dot in the open symbols. A power-law 
fit to these data and its intersection with the green horizontal line in the figure gave the 
parameters listed in Table [31 For 2 < Pr < 4 there also are sufficient data to warrant a 
power-law fit, but there are no classical-state data. We carried out the same analysis as 
for 1 < Pr < 2, and used the result for the classical state for 1 < Pr < 2 as the baseline to 
determine Rat. One sees from the Table that the exponents are, within their statistical 
errors, consistent with 0.38, a value deemed typical of the ultimate state. The values 
of Rat, however, in our view are too low to correspond to the shear- induced boundary- 
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layer transition to turbulence that is expected to be characteristic of the ultimate-state 
transition. Further, as noted in Ref . |T7] , we see that the larger Pr value yielded a lower 
value of Rat, which is opposite to expectations for the BL shear instability. However, 
this trend is not confirmed by the data for the larger Pr range 4 < Pr < 8 (purple open 
triangles), which do not have enough points to warrant an independent power-law fit 
but which are seen to fall between the other two data sets. Thus we must conclude 
that the Grenoble data do not establish an unambiguous trend of Ra^ with Pr. An 
explanation of the different values of Ra* for the different data sets which seems likely 
to us is systematic errors of the equation of state or the transport properties near the 
critical point of helium that were used to calculate Ra and Nu; these errors could well 
change as the pressure and (and thus Pr) are changed. 

The Grenoble data for Pr = 1.1 and 1.3 are of sufficiently high precision and 
sufficiently detailed and plentiful to warrant a closer examination, as is done in Fig. [7] (a). 
There they are compared with our results, which are shown in Fig. [7] (b) on vertical and 
horizontal scales with the same resolution as in (a) but with the horizontal axis shifted 
by two decades. A remarkable difference between the two data sets is that the Grenoble 
data reveal a sharp transition suggestive of a continuous, or supercritical, bifurcation, 
whereas our data show a transition range of about 1.5 decades with indications of multi- 
stability in that range which is inconsistent with a continuous transition. We believe 
that a transition to a turbulent BL is unlikely to be sharp and continuous for at least two 
reasons. First, even for a uniform laminar BL the transition to turbulence does not occur 
at a unique value of the applied stress but rather will depend on the particular prevailing 
perturbations. In the time average this should lead to some rounding of the observed 
transition. Second, in the RB case the laminar BLs are not uniform. Rather, due to 
plume emission, they are highly fluctuating systems. In addition, they are non-uniform 
on longer length scales in the horizontal plane [39] when a LSC is present. Because of 
their spatial inhomogeneity they are unlikely to undergo a simultaneous transition from 
the laminar to the turbulent state at all lateral positions. These inhomogeneities are 
consistent with the existence of a transition range, as observed by us. 

Finally, in Fig. [6l we show as purple open circles with plusses and error bars the 
results obtained from direct numerical simulation by Stevens et al. [80] . These data are 
for Ra < 2 X 10^^. They do not show a transition for Ra < 2 x 10^^ to a state with 
a larger effective exponent, and in that sense differ from the Grenoble data but agree 
with the Oregon data and with our results. 

5.2.2. Comparison with the Chicago data: The "Chicago" data are the earliest 
measurements of Nu at very large Ra. For P = 0.5 they were reported by X. Z. Wu 
[79] . These results, after a correction [HI] which shifted the data upward by an amount 
that varied smoothly from about 6 % near Ra = 10^ to about 10% near Ra = lO^'^, are 
shown in Fig. |8](a) as open black circles. Their trend with Ra suggests 7e// — 0.29 [74] . 
However, they had not been corrected for the finite sample- wall conduction [711 175] . 
That correction, based on Model 2 of [71] , was applied in Ref. [71] and yielded the open 
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Figure 8. Comparison of the present results (solid black circles) with (a) the Chicago 
data |79], (b) data from [14], and (c) to (e) data from [17]. We used different symbols 
for different ranges of Pr. Red open circles: Pr < 1. Green open squares: 1 < Pr < 2. 
Blue open diamonds: 2 < Pr < 4. Purple open triangles: 4 < Pr < 8. The horizontal 
black dashed line is our best estimate for the classical state. The other lines are fits 
described in the text. 
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Table 3. Parameters obtained from power-law fits to various data sets. 
Ra„Mn, Ramax, Pfmrn, and Pimax are the limits of the data sets used in the power- 
law fits. Nu^g^ is the reduced Nusselt number Nu/Ra°'"^^^ used in the classical regime 
to determine the transition Rayleigh number Ra^, which is taken to be the intercept 
with the power-law fit above Rat. 



purple squares in the figure. Those corrected data are consistent with 7e// = 0.312, but 
the pre-factor Nq of a power-law fit is about 10 % lower than indicated by our results. 

Just above Ra = 10^^ there is a discontinuity in the data, but this apparent 
"transition" is not to a state with a larger 7e//- We believe that this phenomenon 
is associated with a change of the place in the phase diagram near the critical point 
of helium where the data were taken, and that it is due to errors in the equation of 
state rather than a genuine change in the dependence of Nu on Ra. Thus we conclude 
that the Chicago data do not reveal any evidence of a transition to the ultimate state 
below their largest Ra ~ 3 x 10^^. Given the increase of Pr with Ra for these data and 
the expected dependence of Ra* on Pr, we find that the absence of a transition to the 
ultimate state in these data is consistent with the Oregon data and with our results. 

5.2.3. Comparison with data from Roche et ai. Finally, in Fig. Mjo) to (e), we display 
several data sets obtained in different sample cells, all with F = 0.5, by Roche and 
co-workers [HI ^7\. Many of them show a transition at Ra* from the classical state 
to a state with a larger 7e//- Whenever there are adequate classical-state data, these 
sets are consistent with the exponent 0.312, but the pre-factor of the corresponding 
power law varied between different sets and generally was somewhat lower than our 
result A^o = 0.104. We analyzed the sets with sufficient data both above and below 
Ra^ as described above for the Grenoble data. The results are given in Table [31 The 
effective exponent varies from 0.371 to 0.404. A reason for this variation is not obvious 
to us. The value of Rai varied significantly as well, from 0.9 x 10^^ to 7 x 10^^, again 
for non-obvious reasons. 

In Fig. [9] we show Ra* as a function or Pr as red solid circles. A trend with Pr is, in 
our view, not firmly established, although the point at the largest Pr suggests a decrease 
with increasing Pr. Such a decrease would be contrary to the theoretical expectation 
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Figure 9. The transition Rayleigh number Rat as a function of Pr. The data are from 
Table [H Blue diamond: HPCF-IIe. Red circles: The remaining points in the table. 
The solid line is a fit of the theoretically expected dependence Rat oc Pr^'^^ (Eq. |6]) to 
the HPCF-IIe point. 



for a shear-driven BL transition to turbulence. The theoretically expected result is 
obtained by assuming that the transition to the ultimate state occurs at a Rayleigh 
number Ra* at which the bulk Reynolds number (and thus also the shear Reynolds 
number Re* oc \/Re*)) attains a certain constant value Re*. Taking Re oc Ra^'^^^^/Pr"'™, 
one finds 



Here the constant coefficient was adjusted so that the relation passes through the point 
Ra* = 1.1 X 10^4 for Pr = 0.8 shown as the blue diamond in the figure, which is the result 
from the present work. The line through that point corresponds to Eq. |6l One sees that 
the red circles do not reproduce the predicted increase of Ra* with Pr, as already noted 
by Roche et al. [17]. However, we do not find convincing evidence in these data that 
larger values of Pr favor the transition to the state with the larger exponent. 

5.3. Open sample 

As discussed in Sec. HI initially (for HPCF-I, Ha, and lib) gaps with an average width 
of about one mm were left between the top and bottom plates and the side wall in order 
to permit the SFg to pass from the Uboot into the sample. These samples are referred 
to as the "open" samples. Although the gap widths were negligible compared to the 
sample dimensions, it was appreciated later that they could significantly modify the 
system, depending on the temperature difference — T(/ between the sample (Tm) and 
the Uboot (Tu). When Tm < Tu {Tm > Tu), then the sample density is larger (smaller) 
than the density of the ffuid in the Uboot. In the presence of gravity this density 
difference will lead to ffow through the gaps, in the upward (downward) direction when 
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Figure 10. The reduced Nusselt number Nu/Ra°''^ as a function of the Rayleigh 
number Ra on a logarithmic scale. The purple solid circles, black solid squares, and 
red solid diamonds are from HPCF-I, HPCF-IIa, and HPCF-IIb respectively, using 
SFg at various pressures. The open red circles are from HPCF-I using N2. All of 
those data were shown before in Ref. [23] • The open black and solid green circles 
are new data from HPCF-IIb using SFg at 12.2 bars (Pr = 0.823) and 18.8 bars 
(Pr = 0.863) respectively. The small black dots are the present results for HPCF- 
Ile (the closed sample). The lines correspond to power laws Nu = iVoRa'''''-'^ with 
No = 0.111, Jeff = 0.308 (dashed black line), Nq = 0.674, jeff = 0.25 (red solid line). 
No = 12.3, 7e// = 0.16 (black dotted line). No = 0.0277, 7e// = 0.36 (black solid line), 
and No = 0.104, 7e// = 0.312 (dotted blue line). The labels U, Li, and L2 identify the 
"upper branch", "lower branch 1", and "lower branch 2". 



Tm — T[/ > (< 0). Such an externally imposed flow can be expected to modify the 
ideal RBC system. 

Previously [7S1 ESI 123] we presented some measurements of Nu for the open system. 
These measurements were made with = 25°C. During those early stages of this 
investigation Tu was not actually measured, but more recent experience suggests that 
it was between 23 and 24°C. Thus the data are for — T[/ > 0. The results are shown 
again in Fig. [10] in the form of Nu/Ra°'^ as a function of Ra. The purple solid circles, 
black solid squares, and red solid diamonds are from HPCF-I, HPCF-IIa, and HPCF- 
IIb respectively and were obtained with SFg. The open red circles were measured using 
Nitrogen in HPCF-I. For comparison the closed-sample data are shown as small black 
dots, with the power-law fit in the classical range (see Fig. [3]) with 7e// = 0.312 shown 
as a dotted blue line. For Ra < 4 x 10^^^ the open-sample data are described well by 
a power law with 7e// = 0.308 (the dashed black line), reasonably consistent with the 
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classical state; however, they are somewhat lower than the closed-sample data. 

With increasing Ra the measurements from HPCF-I (purple solid circles in Fig. [T0|) 
revealed a transition to a new state beyond Rati — 4 x 10^^. Within our resolution Nu 
was continuous at that transition, and the transition was sharp. In those respects the 
transition differed from the ultimate-state transition seen for the closed sample but was 
similar to the transitions at Raj = 0(10^^) found by Chavanne et al. [12] (see Fig. [6] 
and[7](a)) and Roche et al. [17]. However, the present case is very different from the 
previous ones in that 7e// decreased whereas for the prior work it increases. Above Rati 
the data could be described well by a power law with 7e// = 0.25, as shown by the solid 
red line in the figure. We shall refer to this state as the "Li" state or branch. 

In a recent paper [66] Grossmann and Lohse pointed out that, with increasing Ra, 
the transition away from the classical state can take several possible forms. Whereas 
previously [5l |6] the ultimate-state transition was assumed to involve a simultaneous 
shear-induced transition to turbulence in both the viscous and the thermal BL, this 
need not be the actual sequence of events. They proposed three possibilities that may 
prevail when the viscous BL becomes turbulent. One of these is that the thermal BL 
remains laminar {i.e. of the Prandtl-Blasius type in the time average) and that the 
heat transport is background dominated (see Sec. III.B of [66j). For that case they 
derive Nu ~ Ra^''^, with logarithmic corrections which yield 7e// — 0.22 to 0.23. Within 
experimental and theoretical uncertainties this is consistent with our result 7e// — 0.25 
for the Li branch. 

Further measurements, with HPCF-IIa (black solid squares in Fig. [TOl) and HPCF- 
Ilb (red solid diamonds in Fig. [TOj) . revealed the existence of yet another branch which 
we labeled as L2. The transition to this branch occurred at Rat2 — 2.5 x 10^^. The 
new branch co-existed with Li. The precise conditions that determined which branch 
was chosen by the system were not explored in detail, but they involved — Tjj and 
thus the strength of the presumed external current entering or leaving the sample. A 
power-law fit to the L2 data yielded 7e// — 0.16. This result is consistent with the 
theoretical value (see Sec. III. A of [66]) 7 = 1/8 which logarithmic corrections for a 
state with a turbulent viscous BL and a laminar thermal BL, with the heat transport 
dominated by the emission of plumes. For this state the logarithms modify 7 so that 
7e// ~ 0.14, not very different from the experimental value. 

When it was realized that T^—Tu played an important role in the choice of the state 
of the system, new measurements were made with HPCF-IIb and with Tm — Tu < —2 K. 
Those data are shown as open black and solid green circles in Fig. [TOl They fall somewhat 
above the closed-sample data (small black dots), but seem reasonably consistent with a 
perhaps slightly modified transition to the ultimate state. This upper branch is labeled 
"f/". A power-law fit to these data yielded ■jeff = 0.36 and is represented by the solid 
black line in the figure. 
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Figure 11. The reduced Nusselt number Nu/Ra°''^ as a function of the Rayleigh 
number Ra on a logarithmic scale. Small purple symbols: open sample (see Fig. 1101) . 
Small black circles: closed sample, T„i — Tu < — 2 K (see Fig. [5]). The larger red 
symbols are for the half-open sample. Solid squares: Tm — Tjj < —2 K. Open squares: 
T,„ ^Tu>2K. 



5.4- Half-open sample 

As a first attempt to prevent the external currents through the sample due to the 
difference between Tm and Tu, we sealed the bottom plate to the side wall but left the 
gap between the top plate and the side wall so that the gas could still enter the sample. 
We refer to this case, which is HPCF-IIc, as the half-open sample. The results for Nu 
are shown in Fig. [H] as solid red squares for Tm — Tu < —2 K and as open red squares 
for Tm — Tu > 2 K. As might be expected, the data for Tm — Tu < —2 K agree well 
with the results from the closed sample (small solid black dots). This is so because the 
relatively dense gas, which in the open sample escapes through the gap at the bottom 
plate, can not do so in this case. However, the case where the sample gas is less dense 
than the gas in the Uboot (Tm — Tu > 2 K, open squares) seems to be influenced by 
external currents near the top plate and the corresponding data are close to those for 
the open sample (purple solid circles). Indeed they even reveal the transition at Rati to 
the Li state. 
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6. Summary 

In this paper we reported results for the Nusselt number Nu over the range of the 
Rayleigh number Ra from 3 x 10^^ to 10^^. Data were presented for three different 
sample cells, all of cylindrical shape and of aspect ratio F = 0.50. The cells, known 
as the High Pressure Convection Facilities (HPCFs), were located in a pressure vessel 
referred to as the "Uboot" of Gottingen. The Uboot and the HPCF were filled with 
the gas sulfur hexafluoride (SFg) at various pressures up to 19 bars. This fluid had a 
Prandtl number Pr which varied over the narrow range from 0.79 to 0.86 as the pressure 
(and thus Ra) changed from its smallest to its largest value. 

One of the samples was completely sealed. In that case a 2.5 cm diameter tube 
penetrated the side wall at mid height and allowed the SFg to enter the HPCF from the 
Uboot; that tube was then closed off by a remotely operated valve after an equilibration 
time of several hours and before measurements were made. This version of the HPCF 
is known as the "closed" sample. Another version was allowed to have a small gap, of 
width approximately equal to one mm, between the side wall and the top and bottom 
plate to allow the sample gas to enter the HPCF. This version is known as the "open" 
sample. A third version, known as the "half-open" sample, had the side wall sealed to 
the bottom plate while the gap was allowed to persist at the top plate. 

It turned out that the three samples produced qualitatively different dependences 
of Nu on Ra. Only the closed sample can be regarded as corresponding well to the 
idealized closed-system Rayleigh-Benard problem. For the open and half-open samples 
gas currents were able to enter and leave the sample through the gaps at the plates, 
driven by the small density differences that existed between the sample fluid at a mean 
temperature Tm and the fluid in the Uboot at a temperature Tu. The results for Nu(Ra) 
were then qualitatively different depending on whether was larger or smaller than T^/. 
Nonetheless the open and half-open samples produced interesting effects. Under certain 
conditions there was a sharp but continuous transition from a state below Ra^ ~ 4 x 10^^ 
where Nu ~ Ra'^'^-'^-'' with 7e// ~ 0.31 to a state above it with jeff — 0.25. Another state 
with 7e// ~ 0.16 was found as well. Possible explanations of these findings in terms 
of different laminar or turbulent conditions in the thermal and viscous boundary layers 
adjacent to the top and bottom plates were offered by Grossmann and Lohse [66]. We 
refer the reader to Sects. I5.3l and [5^ as well as to Ref. |66j for a more detailed discussion 
of these interesting phenomena. 

The primary focus of this paper has been on the closed sample. Even for this case 
the results for Nu(Ra) depended somewhat on Tm — Tu, but in contradistinction to the 
open and half-open samples the dependence was simply a shift of the curve without any 
change in shape. This dependence persisted in spite of the extensive shielding of the 
sample that was provided, and we do not know its origin. However, for Ra < 10^^ we 
found that the effective exponent of Nu was 0.312 ± 0.002, regardless of Tm and T^/. 
This exponent value is consistent with many other measurements at smaller Ra. It also 
agrees quite well with the value 0.323 with follows from a power-law fit to data generated 
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from the Grossmann-Lohse model [M] for Pr ~ 0.8 (albeit for F = 1.00) in the range 
10^^ < Ra < 10^^. Thus we believe that the state of the system with Ra < 10^^ is 
the classical state of RBC with laminar boundary layers below the top and above the 
bottom plate. 

In the range Raj < Ra < Rag, with Raj ~ 1.5 x 10^^ and Rag ~ 5 x lO^'', the system 
gradually underwent a transition and 7e// increased from 0.32 to 0.37 as Ra changed by 
about one and a half decades. We believe that this phenomenon reflects the transition to 
the ultimate state predicted by Kraichnan [5] and re-examined recently by Grossmann 
and Lohse [66]. In the ultimate state, which we found for Ra > Rag, the effective 
exponent was 0.37 ± 0.01. This is reasonably consistent with the predicted asymptotic 
exponent 7 = 1/2 and the expected logarithmic corrections [66]. Our conclusion that 
the state above Rag is the ultimate state is supported by transitions and exponent- 
values found in simultaneous measurements of the Reynolds number Re ^26j (a detailed 
discussion of those results is beyond the scope of this paper). Evidence for the transition 
range, with about the same values of Rajg, can be found also in recent measurements 
of logarithmic vertical temperature profiles that extend throughout most of the sample 

m- 

It is worth emphasizing that the transition from the classical to the ultimate state 
is not a continuous transition where Nu(Ra) is continuous but its effective exponent 
changes discontinuously such as was found at Raj ~ 4 x 10^^ in the open sample. 
Instead there was a transition range spanning a factor of 30 or so in Ra over which the 
transition occurred. In this transition range the results for Nu(Ra) scattered more than 
below or above it and were often irreproducible at the level of our resolution from one 
point to another, suggesting the existence of many states which differed in detail. 

The observed transition range (as opposed to a unique value of Ra*) is not surprising 
for two reasons. First, transitions involving shear-flow instabilities are known to depend 
on the size of prevailing local perturbations and thus will occur at different values of 
the Re (or in our case Ra) number. In addition, in the present case the basic state, i.e. 
classical RBC, is known to have laminar boundary layers which are non-uniform in the 
horizontal plane |39|. Since the BLs and the shear applied to them by the turbulent bulk 
are known to be spatially inhomogeneous, one should not expect a sharp transition. 

In this paper we also provided a re-examination of earlier measurements. Some 
of these indicated a transitions to a state where 7e// assumes a value near 0.38, but 
the transitions occurred at the rather low values of Ra^ between 10^^ and 10^^. All of 
these data involved a range of Pr. We separated them into subsets corresponding to 
a unique value, or spanning only a small range, of Pr. Within each subset one sees 
that the scatter of these data is much less than it is when all Pr values are considered 
jointly. Some of the subsets reveal a well defined transition. These transitions differ 
qualitatively from the one observed by us in that they are sharp and continuous, without 
any resolved transition range. Estimates of the shear Reynolds number for them are in 
the range near 100 or less, which to us seems too low for a turbulent boundary-layer 
shear transition. Based on our measurements of Re, and on recent calculations by DNS 
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of the shear Reynolds number Re^ ^2] , as well as on earlier estimates [2Z] , indicate that 
the boundary-layer shear transition at Ra* to the ultimate state should not occur until 
Ra = 0(10^"^) is reached. Thus, in our view, the transitions observed near Ra = 10^^ to 
10^^ are unrelated to the ultimate-state transition; but we can not offer an alternative 
explanation for their existence. 

This paper concludes with a number of Appendices which examine the possible 
effect of several experimental factors on our results. None of them are found to have a 
significant influence. 
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Appendix A. Effect of tilting the sample 

All data for this study were taken with HPCF-IIe (the closed sample). Measurements 
were made with a leveled sample (a tilt angle /3 < 0.1 mrad, run 100918 to 110102), 
and with the same sample tilted relative to gravity. The tilt was used to localize the 
azimuthal orientation of the large-scale circulation plane. Based on measurements with 
smaller samples of aspect ratio F = 1.00 [32] and 0.50 [83j, both with Pr ~ 4, we do not 
expect a measurable influence of the tilt on the heat transport for our case of F = 0.50 
and Pr ~ 0.8. This is indeed born out by the data. 

In one case the tilt angle was /3 = 3 mrad (run 110115 to 110610). The tilt 
direction was at an angle Op = ?)'k/2 rad relative to the azimuthal origin of the sidewall 
thermometers that were used to monitor any large-scale circulation. In the other case 
we had /3 = 14 mrad and Of^ =? rad (run 110618 to 110919). The results for the Nusselt 
number for the data with — Tu< — 2K are shown in Fig. lAll as blue circles (no tilt), 
green diamonds (/3 = 3 mrad) and red squares (/3 = 14 mrad). There is no signiflcant 
difference between the three data sets, showing that the effect of the tilt on Nu is well 
below one percent. 

Appendix B. Non-Boussinesq effects 

Non-Oberbeck-Boussinesq (NOB) effects [HH [34] on Nu in a gas have been studied 
quantitatively using ethane gas as the fluid [HH IHS]. Although variations of all fluid 
properties contribute in principle, for gases the primary contribution comes from the 
temperature dependence of the density and can be described approximately by the 
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Figure Al. The reduced Nusselt number Nu/Ra^ '^^'^ as a function of the Rayleigh 
number Ra for Tm — Tjj < — 2K as a function of the Rayleigh number for three different 
tilt angles /3 of the sample relative to gravity. Blue circles: "level" , f3 < 10"^ rad. Green 
diamonds: f3 ~ 0.003 rad. Red squares: (3 ~ 0.014 rad. 



parameter a AT (see Fig. 4 of [5^). AH our measurements were made with AT < 21K 
where we expect NOB effects to be smalL For a perfect gas one has aAT = AT/T 
with T in K, which for AT ~ 20K and at a mean temperature of 21°C or 294 K is 
close to 0.07. However, the properties of SFg in our pressure range up to 19 bars show 
large deviations from those of a perfect gas. Figure IBlb shows aAT as a function of 
Ra for the data with T^ — Tjj < —2 K. The different-colored circles are, from left to 
right, for different isobars with pressures of approximately 4, 5, 8, 12, 15, and 19 bars. 
For the points shown as red squares we had 19 < AT < 21 K, which were the largest 
temperature differences employed in this work. The dotted line is the approximate 
perfect gas result. One sees that aAT was mostly below 0.1 and always below 0.14, but 
at the highest AT and pressure exceeded the perfect gas value by about a factor of two. 

In Fig. IBlb we show the reduced Nusselt number for the data with Tm — Tu < —2 K. 
All data are shown as solid black circles, except for the ones with 19 < AT < 21 K which 
are given as red squares. One sees that the points with large AT are systematically lower 
than the remaining results, but only by about 1.4 percent at the highest pressure where 
a AT is largest and by about 0.7 percent at the lowest pressure where a AT is about a 
factor of two smaller. It is surprising that NOB effects reduce Nu in the present case 
because in the case of ethane [HH [HH] they enhanced Nu, albeit by only a small amount. 
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Figure Bl. a): The parameter a AT as a function of Ra. Tlie circles of various 
colors are, from left to right, for approximately 4, 5, 8, 12, 15, and 19 bars. The red 
squares are points for which 19 < AT < 21 K. The dotted horizontal line indicates the 
perfect-gas value 0.071 for AT = 21.0 K and T„ = 294.2 K. b): The reduced Nusselt 
number for T™ — Tjj < —2 K as a function of the Rayleigh number. Black circles: 
AT < 19K. The red squares are points for which 19 < AT < 21K and correspond to 
the red squares in a). 



A reduction of Nu due to NOB effects was observed, however, in the case of a hquid, 
namely water [86], but again only by a percent or so for AT ~ 20 K. In any case, what 
matters for the present work is that NOB effects are negligible for nearly all of our data. 

Appendix C. The parameter ^ of Niemela and Sreenivasan 

Recently Niemela and Sreenivasan [35] introduced the parameter 

where A is the fluid conductivity, r] is the shear viscosity, and A indicates the difference 
of the value of the property following it at the bottom and at the top of the sample. 
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Figure CI. The parameter ^ introduced by Niemela and Sreenivasan [35] (see text). 
Tire data are for the following approximate pressures: blue diamonds: 4 bars; red 
diamonds: 5 bars; green circles: 8 bars; green triangles: 10 bars; red squares: 12 bars; 
red triangles: 14 bars; yellow circles: 15 bars; green squares: 16 bars; black circles: 18 
and 19 bars. 



Thus, is a measure of the relative sizes of the contributions to NOB conditions from 
A and rj. The authors state that their data for a cell with T = 1.0 [191 ES] and T = 4.0 
[20] make a transition from a state of lower to one of higher Nu as ^ changes in their 
experiments from positive values of order one to negative values of -1 or less, both 
states having a common scaling exponent with a value slightly above 0.3. The authors 
assert that the transition between the two states, where roughly Nu ~ Ra'''^, should not 
be interpreted as the ultimate state; rather, the data correlate well with the material 
parameter ^. 

The values of ^ for our data are plotted in Fig. ICll We see that at a given pressure 
^ is essentially independent of Ra, and thus of AT. The value of ^ depends slightly 
upon the pressure, varying over the range 1.6 > ,^ > 1.3 as the pressure changes from 
about 4 to about 19 bars. Thus, for all of our data ^ remains positive and slightly above 
one and the transition seen in the data of [T9| 120] [35] , which were taken near the critical 
point of helium, is not expected to occur in our fluid-property range . 

Appendix D. Effect of horizontal temperature variations in the top and 
bottom plate 

Maintaining the top and bottom plate at a uniform temperature, even in the presence of 
the large vertical heat currents up to a kW or so, is a significant experimental challenge, 
especially when the plates become very large as in our experiment where the diameter 
was well over a meter and the mass of the plates was of order 200 kg. 

Our top plate was cooled by water passing through a pair of double spiral channels 
[23] . The two members of the pair were in parallel, and each pair consisted of anti- 
parallel {i.e. an inward and an outward flowing) spirals. The distance between the 
centers of adjacent channels was 2.54 cm, and the channels had a depth d = 1.26 and a 




Figure Dl. a.) The normalized radial temperature variation (T^ — To)/AT as a 
function of the temperature difference Tm — Tjj between the mean sample temperature 
and the Uboot temperature (for details, see text). Red symbols: top plate data for 
the weaker water circulation driven by a Nelab RTE7 circulator. Blue symbols: Data 
for the top plate with the stronger water circulation driven by a Wilo pump. Black 
symbols: Data for the top member of the bottom-plate composite. Open circles, 
solid circles, open squares, and solid squares are for Tq — Ti with i = 1, 2, 3, and 4 
respectively, b.) The azimuthal temperature variation for two examples. The solid 
lines are a fit of a sine function with an adjustable phase and amplitude to the data. 



width w = 1.27 cm (cross sectional area of 1.60 cm^). These dimensions were kept small 
to minimize lateral temperature variations on the small length scales of order the plate 
thickness. However, as will be seen below, the small channel cross section did lead to a 
significant flow resistance, and a somewhat larger channel cross section might have been 
optimal. There is a constraint on the channel depth provided by the plate thickness, 
and a much greater thickness would lead to excessive weight and cost of the plate. 

The plate temperatures were determined by five thermometers in each of the top 
plate (between the water channels) and the top of the bottom-plate composite [53] . All 
ten thermometers were placed in small holes in the plates, with their tips within a mm or 
so of the copper-fluid interface. One thermometer (Tq) was located at the plate center. 
The other four Tj, i = 1, ...,4, were equally spaced azimuthally at angular intervals of 
Tc/2 rad and were positioned radially at a distance of 0.42D from the center. The plate 
temperatures and Tf used to calculate AT (and thus Nu and Ra) were the averages 
of the flve readings in a given plate. 

Here we report on results obtained with HPCF-IIe. Initially (run 110918 to 110410) 
the water flow was generated directly by a Neslab RTE7 circulator capable of generating 
a pump head corresponding to approximately 0.5 bars which yielded a relatively small 
flow rate near 15 cm^/s (0.05 m^/h) in each of the two double spirals, corresponding to 
a mean velocity v of about 10 cm/s or a Reynolds number Re = vd/ v = C(IO^). Judged 
by the value of Re, the flow probably was laminar. In Fig. IDlb we show (Tj — Tq)/AT, 
z = 1, 4 as a function of Tm — Tjj as red symbols. Those data were taken at a pressure 
of 18.8 bars and with AT ~ lOK (run 110311 to 110326) corresponding to Ra ~ 8 x W^, 
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Figure D2. The reduced Nusselt number as a function of tlie Rayleigh. number for 
Tm — Tjj < — 2K. Open red (solid blue) circles: data taken with low (high) water flow 
rate in the top-plate cooling channel. 



but similar measurements over a range of AT and at other pressures gave similar results. 
The dependence on — Tu is weak. The azimuthally averaged radial temperature 
variation ((Tj)j=i_..._4 — To)/AT is a little over two percent, which was judged larger than 
desirable. The azimuthal variation of (Tj — Tq)/AT for — Tu = — 3.2K (run 110311) 
is shown in Fig. IDlb as red symbols. It can be fit well by a sine curve, suggesting that 
its origin is the large-scale circulation in the bulk fluid [291 EH 188] . 

In order to increase the flow rate of the water in the top plate, we modified the top- 
plate water cooling circuit (starting with run 110426). We drove a primary circuit by a 
Wilo model MHI 205-1 pump which could generate a pressure differential across the top- 
plate water channel of 4.2 bars which yielded a flow rate 170 cm^/s (0.6 m^/h) in each 
double spiral, corresponding to a mean flow velocity v ~ 100 cm/s or Re = 0{1Q^). This 
flow almost certainly was turbulent. The primary water circuit was thermally coupled 
by a heat exchanger to a secondary cooling circuit driven and temperature-controlled 
by the Neslab RTE7 circulator. This circulator was servoed so as to maintain Tt at its 
desired value. The results for (T— To) / AT are shown in Fig. IDlk as blue symbols and for 
two values of T^ — Tjj. There was a signiflcant improvement, with ((Tj)j=i^... 4 — To)/AT 
as small as 0.8 percent. The remaining azimuthal variation for T^ — Tu = —2.2 K is 
shown as blue symbols in Fig. IDlb . 

In order to judge whether one or the other of the top-plate cooling circuits is 
adequate to yield valid results for Nu, we show in Fig. ID2I results for Nu(Ra) obtained 
by the flrst method as red open circles. Data obtained with the improved top-plate 
cooling method are represented by the solid blue circles. The two sets agree with each 
other very well and essentially within their scatter, indicating that there is no systematic 
dependence on the cooling method and that either method was indeed adequate for 
Nusselt number measurements. 
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The bottom plate was a composite of two copper plates with a thin Lexan layer 
between them. The heat current passing through this composite was generated by Joule 
heating with a heater wire imbedded in grooves milled into the underside of the bottom 
member of the composite and was uniformly distributed over the entire sample area 
[23] . The current had to pass through the Lexan thermal barrier before passing through 
the top copper member and then entering the fluid. Results for (Tj — Tq)/AT obtained 
for runs 110918 to 110410 are shown as black symbols in Fig. ID lb . Their average values 
are close to zero, showing that the radial temperature variation in the bottom plate was 
remarkably small and indeed negligible. 

Appendix E. Evidence for a closed sample 

HPCF-IIe was supposed to be completely sealed, except for a 2.5 cm diameter pipe 
entering it at half-height through the side wall. This pipe terminated outside the sample 
in a remotely operable valve which was to be open during an equilibration period to 
allow pressure equilibration between the Uboot and the sample, but closed when the 
system had reached a steady state after a new set point of the top- and bottom-plate 
temperatures had been established. In Fig. lElk we show the pressure inside the sample 
during this equilibration period. With the valve open, the pressure of the Uboot and 
sample still drifted at a significant rate. In this case, however, the valve was closed 
prematurely. One can see that the sample pressure rose. When the valve was opened 
again, the pressure returned back to the Uboot pressure which all along had been drifting 
downward, albeit at a lesser rate. We regard this measurement as evidence that sealing 
the sample was indeed successful. 

Figure lElb shows the pressure with higher resolution after a steady state had been 
reached. With the valve open, the noise level of the readings, reflecting the instrumental 
noise, is quite small but noticeable. After closing the valve, the noise level increases 
because the genuine noise in the pressure due to the turbulent convection in the sample 
becomes noticeable because it can no longer be smoothed out by exchange of fluid with 
the Uboot. This phenomenon also indicates that sealing the sample was successful. 

Appendix F. Effect of side-shield mismatch with the sample temperature 

As shown in Fig. [1], the side wall of the sample is surrounded by a side shield ("SS" in 
the flgure) which is intended to prevent heat loss or input through the side wall. The 
temperature of this shield was controlled at the mean temperature T^. Since the sample 
center temperature will differ slightly from due to non-Boussinesq effects, it is 
important to ask whether heat input to the side wall due to the temperature difference 
Tm — Tc is signiflcant. Further, in view of the large size and mass of the shield, it is 
appropriate to ask how successful its temperature control was at a specifled setpoint. 
Thus, a run was conducted in which the shield temperature was deliberately servoed 
at a displacement from which is much larger than — or any deviations from 
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Figure El. a): The sample pressure as a function of time. The data were taken on 
Sep. 14 2010, during a time interval when the Uboot temperature and pressure had 
not yet fully equilibrated. The sharp increase at the time indicated by the first arrow 
occurred when the valve was closed, and the sudden decrease at the second arrow was 
caused by opening the valve again, b): The sample pressure as a function of time since 
the beginning of run 100918. Initially, the valve was open and the sample pressure could 
equilibrate with the Uboot pressure. When the valve was closed a.t t = 40200 s (the 
location of the arrow) , the pressure could no longer equilibrate with the Uboot pressure 
and the fluctuations due to the turbulent convection became noticeable. This run is 
for T„ = 20.69°C, 18.4 bars, and AT = 11.37K, corresponding to Ra = 9.25 x 10^'^. 




Figure Fl. The fluid temperature at a radial distance of one cm from the side wall 
as a function of the height z/L above the bottom of the sample. The red circles (black 
squares) are for (Tss) = {{Tss) = T™ + 1 K). 
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Table Fl. Effect of a side-shield temperature-deviation from tlie mean temperature 
for run 120415. (Tss) is tlie set-point of the shield temperature. Tss{iL/A), i = 1, 2, 3, 
are the measured shield temperatures at the three vertical positions iL/4 as measured 
from the sample bottom. 

Quantity {Tss) = Tm {Tss) =T„, + IK 



T 


21.587 


21.599 


Tc 


21.398 


21.441 


AT 


10.163 


10.191 


Tss{L/A) 


21.629 


22.621 


Tss{2L/4) 


21.602 


22.602 


Tssi^L/A) 


21.550 


22.571 


Q [Watts] 


320.6 


319.6 


Ra 


6.70 X 10^^ 


6.76 X 10 


Nu 


4779 


4747 



perfect temperature regulation, namely at T^ + LO K. Various measured quantities with 
both setpoints are shown in Table [FTI This example is for a nominal AT = 10.1 K. 

First we call attention to the shield temperatures Tss at the three vertical positions 
L/4, L/2, 3L/4 (the sample height L was 2.2 m) measured from the bottom plate. They 
span a vertical distance of over one meter. One sees that the temperature gradient in 
the shield was typically about 0.05 K/m, which suggests that the shield temperature 
was uniform throughout within about 0.1 K. 

A small shift of Tm and AT can be seen to have occurred due to the shift of the 
shield temperature by one K. This could happen because the bottom plate is a composite 
with a lower (BPb) and an upper (BPt) member (see Fig. [1]), with the temperature of 
BPb controlled. Thus the temperature of the BPt could adjust itself to prevailing heat 
currents and influence T^ as well as AT. The effect of a one K change in the shield 
temperature caused a change of 28 mK, or 0.3%, of AT and a shift of Tm by 12 mK. The 
effect of imperfect shield regulation on these quantities can be regarded as negligible. 

The Rayleigh number was shifted by about 0.9%. This is due in part to the change 
of AT and in part due to a change of the fluid properties associated with the shift of 
Tm- Since actual temperature offsets are much less than one K, this effect is not serious. 
Similarly Nu is affected, by about 0.67%. This shift is due both to the change of AT 
and to a small change of the heat current needed to maintain the temperature of the 
bottom member of the bottom-plate composite (BPb) at the specified set-point. Again 
this is not a serious shift. 

For completeness we also discuss the influence of the shield temperature on the 
sample temperature near the side wall. Although not directly relevant to the present 
paper, it is important for other related investigations of temperature profiles in the bulk 
of the sample ^3] . For that purpose eight thermometers were mounted in the sample at 
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a radial position one cm from the side wall and at eight vertical positions Zj. The eight 
temperatures Q{zj) = [{T{zj)) —Tm]/(Th — Tt) (we denote the time average by (...)) are 
shown in Fig lFll as a function of z/L. One sees that there is no change in the shape of 
the vertical temperature profile, but there is a small shift as had been indicated already 
by the shift of Tc shown in the table. 

Finally we note that there are several other thermal shields in the system (see 
Sec. H]), but we believe that their influence on the system performance is smaller than 
that of the side shield. 



Heat transport by turbulent convection for Ra < 10 
Appendix G. Data tables 



37 



Table Gl. SFg, HPCF-IIe, level sample (tilt angle /3 < 0.0001 rad). 



Run No. 


P (bars) 


T 


(°C) 


Tu' 


(°C) 


AT (K) 


Ra 




Pr 


Nu 


100918 


18.419 


20, 


.691 


23. 


,777 


11.370 


9.25IIe+14 


0, 


,861 


5507, 


.45 


100924 


18.466 


20, 


.788 


24. 


,286 


11.564 


9.477C+14 


0, 


,862 


5579, 


.63 


100925 


18.502 


21, 


.844 


24. 


,446 


11.679 


9.167e+14 


0, 


,862 


5447, 


.81 


101004 


8.011 


20, 


.896 


24. 


,343 


11.786 


6.086e+13 


0. 


,799 


2145, 


.55 


101025 


8.000 


20, 


.973 


23. 


,993 


2.947 


1.514e+13 


0. 


,799 


1371 


.90 


101027 


8.003 


20, 


.950 


24. 


,101 


4.896 


2.519e+13 


0. 


,799 


1616, 


.40 


101029 


7.998 


20, 


.925 


23. 


,947 


7.846 


4.03IIe+13 


0. 


,799 


1879, 


.88 


101031 


8.005 


20, 


.895 


24. 


,160 


11.786 


6.075e+13 


0. 


,799 


2146, 


.16 


101102 


8.016 


20, 


.870 


24. 


,549 


17.732 


9.179e+13 


0. 


,799 


2449, 


.66 


101114 


11.751 


20, 


.977 


24. 


,611 


13.447 


2.063e+14 


0, 


,820 


3224, 


.14 


101116 


11.762 


21, 


.021 


24. 


,849 


20.030 


3.079e+14 


0, 


,821 


3675, 


.61 


101108 


14.891 


20, 


.920 


25, 


,015 


17.826 


6.060C+14 


0, 


,842 


4659, 


.07 


101109 


14.894 


21, 


.247 


25, 


,008 


18.479 


6.21IIe+14 


0. 


,842 


4694, 


.89 


101120 


14.894 


21, 


.382 


24. 


,994 


18.749 


6.27IIe+14 


0. 


,842 


4710, 


.51 


101122 


14.843 


21, 


.498 


24. 


,246 


11.988 


3.945e+14 


0. 


,841 


4045, 


.17 


101124 


14.810 


21, 


.272 


23. 


,786 


7.938 


2.61IIe+14 


0. 


,841 


3526, 


.58 


101126 


14.792 


21, 


.282 


23. 


,526 


4.960 


1.624e+14 


0. 


,841 


2995, 


.86 


101130 


14.802 


20, 


.970 


23, 


,740 


9.934 


3.298C+14 


0, 


,841 


3814, 


.08 


101202 


14.839 


21, 


.011 


24, 


,271 


15.011 


5.021C+14 


0, 


,841 


4382, 


.99 


101204 


14.898 


21, 


.517 


25, 


,061 


19.021 


6.339C+14 


0, 


,842 


4724, 


.05 


101207 


14.873 


20, 


.978 


24, 


,757 


17.942 


6.059C+14 


0, 


,842 


4661, 


.48 


101214 


18.739 


20, 


.933 


23. 


,730 


11.523 


1.004e+15 


0. 


,863 


6154, 


.13 


101215 


18.769 


20, 


.971 


23. 


,918 


9.934 


8.699e+14 


0. 


,863 


5410, 


.40 


101217 


18.755 


20, 


.973 


23. 


,781 


7.945 


6.93IIe+14 


0. 


,863 


4975, 


.08 


101218 


18.814 


21, 


.644 


24. 


,306 


12.274 


1.051e+15 


0. 


,863 


5871, 


.83 


101219 


18.817 


21, 


.612 


24. 


,295 


11.214 


9.620C+14 


0, 


,863 


5599, 


.07 


101221 


18.812 


21, 


.536 


24, 


,246 


10.065 


8.657C+14 


0, 


,863 


5375, 


.00 


101223 


18.785 


21, 


.566 


23, 


,998 


9.125 


7.785C+14 


0, 


,863 


5183, 


.77 


101225 


18.781 


21, 


.583 


23, 


,960 


8.160 


6.949C+14 


0, 


,863 


4981, 


.14 


101227 


18.751 


21, 


.490 


23. 


,682 


6.975 


5.924e+14 


0. 


,863 


4710, 


.15 


101229 


18.709 


21 


.432 


23. 


,289 


0.861 


7.258e+13 


0. 


,863 


1754, 


.60 


101231 


18.740 


21, 


.483 


23. 


,570 


5.961 


5.051e+14 


0. 


,863 


4448, 


.59 


110102 


18.735 


21, 


.487 


23. 


,522 


4.971 


4.205e+14 


0. 


,863 


4160, 


.20 
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Table G2. SFe, HPCF-IIe, First tilted sample, tilt angle /3 = 0.003 rad. 



Run No 


P (l)arsi 


m 


(°C) 


u 




AT (K) 


Ra 




Pr 


Nu 


110115 


4.181 


20. 


,950 


23. 


,734 


4.896 


5 I76e+12 


0. 


,787 


966.79 


110117 


4.187 


20. 


,931 


24. 


,142 


19.859 


I07e+13 


0. 


,787 


1497.77 


110129 


4.180 


21. 


,083 


23. 


,621 


3.164 


3 334e+12 


0. 


,787 


843.02 


110131 


4.187 


20. 


,930 


24. 


,159 


19.858 


2 107e+13 

^ • X \j f 1 X kj 


0. 


,787 


1498 53 


110202 


4.184 


20. 


,855 


23. 


,987 


13.708 


1 455e+13 


0. 


,787 


1335 73 

X KJKJKJ • t KJ 


110204 


4.184 


20. 


,895 


23, 


,947 


9.791 


1 038e+13 


0. 


,787 


1201.21 


110206 


4.185 


20. 


,938 


24. 


,046 


5.875 


6 228c+12 


0, 


,787 


1024.77 


110221 


1 8 565 


21. 


040 


24. 


198 


12 067 

X ^ . \J\J 1 


1 nOlr-l-15 


n 


862 


5707 1 6 

kJ t \J t . X \J 


110227 


1 8 574 


21. 


034 


24. 


286 


1 2 057 


1 nnTTp_|_i5 


n 


862 


571 6 32 

tJ t J-\J • 


110311 


18 574 


21. 


034 


24. 


260 


1 059 







862 


5302 69 


1 10417 


1 8 603 


21. 


036 


24. 


563 


7 070 






KJ t 


862 


4694 05 


110418 


18.572 


21. 


,035 


24. 


,259 


5.068 


4 211e+14 


0. 


,862 


4133 38 


1 1 H/l 1 Q 


io.004 


on 
zu. 


QQfi 

,yyD 


0/1 

Z'i. 


,uyo 


o.yyu 


o.ou / e-|-i4 


u. 


SRO 
,oDZ 


Q7i^S 70 
O / Oo. / Z 


110420 


18.573 


20. 


,994 


24, 


,283 


5.990 


4.990e+14 


0. 


,862 


4412.69 


110421 


18.620 


22, 


,003 


24, 


,600 


6.004 


4.813e+14 


0, 


,863 


4374.31 


110524 


11.837 


24. 


,996 


28, 


,197 


13.985 


1.956e+14 


0. 


,820 


3157.24 


110525 


11.827 


23. 


,264 


28, 


,250 


15.655 


2.294e+14 


0, 


,821 


3336.36 


110526 


11.828 


23. 


,974 


28. 


,184 


15.941 


2.289e+14 


0. 


,821 


3331.77 


110526 


11.820 


23. 


,959 


28. 


,025 


11.912 


1.708e+14 


0. 


,820 


3023.89 


110527 


11.809 


23. 


,961 


27. 


,803 


7.918 


1.13IIe+14 


0. 


,820 


2633.26 
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Table G3. SFe, HPCF-IIe, Second tilted sample, tilt angle /3 = 0.014 rad. 



Run No 




^ m 




u 




AT (K) 


Ra 




Pr 


Nu 


110625 


4.971 


24, 


.017 


27. 


,841 


1.033 


1 533e+12 


0. 


,789 


661.28 


110629 


6.953 


23, 


.855 


28. 


,251 


19.704 


6 607e+13 

W.WWI \j \ X (J 


0. 


,795 


2178, 


.61 


110630 


6.946 


23, 


.908 


27. 


,964 


9.814 


3 279e+13 


0. 


,795 


1744, 


.85 


110701 


6.943 


23, 


.953 


27. 


,841 


4.904 


1 635e+13 

X .W(JWvy 1 X kJ 


0. 


,795 


1395, 


.79 


110716 


18.237 


20, 


.992 


23. 


,753 


4.981 


3 817e+14 


0. 


,861 


4011, 


.57 


110717 


18.275 


20, 


.998 


24. 


,135 


9.990 


7 725c+14 

i t i ^'.Jkj \ XT: 


0. 


,861 


5162, 


.00 


110718 


18.282 


20, 


.997 


24. 


,244 


12.980 


1 006c+15 

x.www*^ 1 xw 


0. 


,861 


5707, 


.97 


1 10718 


IS 303 


21, 


.524 


24. 


363 


1 3 033 

X O . WOO 


q SOAp 1 14 


n 

w . 


861 

1 w W X 


5629, 


.48 


110720 


1 8 383 


21, 


531 

. '~J X 


25 


,210 


1 3 647 


X . w wwv^n^ J- 


n 

w . 


862 


5790, 


.57 


110722 


18 379 


21 


530 


25 

^ W • 


156 

. X(J w 


13 046 

X (J. W^W 


1 n09p+15 

X . w w «_/ v^n^ J- w 




w. 


862 


5673, 


.03 


1 1 0723 


1 8 379 


21, 


531 

. W(J X 


25 

^ W • 


135 


1 2 050 

X ^ . WWW 


(_/ . (J X W 1 X j: 




w. 


862 

. UW^ 


5496, 


.51 


110724 


18.384 


21, 


.523 


25. 


,183 


11.037 


8 548e+14 


0. 


,862 


5323, 


.81 


1 1 0725 

X X\J 1 


18 381 

X U. (JU X 


21, 


515 

. W X (J 


25 

^ w • 


,144 


1 023 

X w. W^ w 


f • f W kJ\J \ X t: 




w. 


862 

. UW^ 


5144, 


.72 


110726 


18.370 


21, 


.512 


25. 


,037 


9.019 


6 q65e+14 

W . i_/ W W ^ 1 X ^ 


0. 


,861 


4950, 


.18 


110727 

X X W 1 ^ 1 


1 8 357 


21, 


507 


24, 


906 


8 007 

O . WW 1 


6 166r+14 

W . X WW^ 1 X ^ 


n 

w . 


861 


4741, 


.89 


1 10728 

X X W 1 


1 8 332 


21, 


500 


24. 


654 


6 000 

w . www 


4 59TTr+14 

t: . • .7 ly X XLy 1 X t: 


n 

w ■ 


861 

1 W X 


4264, 


.08 


1 10731 

X X w 1 ox 


1 8 340 


21, 


498 


24. 


760 

1 1 \J\J 


4 993 


3 830P+14 

W .wWW^ 1 X ^ 


n 

w . 


861 

1 W X 


3973, 


.54 


110801 


18.323 


21, 


.004 


24. 


,672 


5.004 


3 qi6e+14 

KJ • iJ X W V. 1 X ^ 


0. 


,861 


3980, 


.43 


110802 


18.341 


21, 


.001 


24. 


,848 


6.997 


5 501e+14 

W • W W X \j 1 X j: 


0. 


,861 


4531, 


.63 


110803 


18.339 


21, 


.261 


24. 


,791 


5.520 


4 281e+14 


0. 


,861 


4140, 


.67 


1 1 0803 


16 166 

X W. X WW 


21, 


475 


24. 


995 


5 945 






w. 


850 

. uw w 


3521, 


.90 


110804 


16.168 


21 


.475 


25. 


,014 


5.946 


2 697e+14 

^•Wt^ f V? 1 X^ 


0. 


,850 


3521 


.17 


110805 


16.173 


21, 


.467 


25. 


,078 


8.431 


3 831e+14 


0. 


,850 


3996 


.65 


110806 


16.182 


21, 


.473 


25. 


,197 


11.937 


5 435e+14 

W.^OW^ 1 XT: 


0. 


,850 


4512, 


01 

. W X 


110807 


14.105 


21, 


.649 


24. 


,915 


6.295 


1 720c+14 


0. 


,836 


3035, 


.56 


i iUoUo 




91 


.oou 


94 


QDQ 


P. 9Q1 


1 71Qr.-l-1/l 


n 
u. 


,000 


3039, 


.94 


110809 


14.113 


21, 


.742 


25. 


,037 


8.982 


2.45IIe+14 


0. 


,836 


3431, 


.72 


110810 


14.136 


21 


.899 


25. 


,354 


12.790 


3.49IIe+14 


0. 


,836 


3863, 


.76 


110811 


12.061 


21, 


.453 


25. 


,045 


5.906 


9.694e+13 


0. 


,822 


2503, 


.66 


110812 


12.060 


21, 


.454 


25. 


,034 


5.908 


9.695e+13 


0. 


,822 


2500, 


.66 


110813 


12.065 


21, 


.442 


25. 


,135 


8.380 


1.378e+14 


0. 


,822 


2822, 


.21 


110814 


12.067 


21, 


.434 


25. 


,177 


11.861 


1.951C+14 


0. 


,822 


3171, 


.61 


110815 


12.078 


21, 


.445 


25. 


,383 


16.879 


2.784e+14 


0. 


,823 


3564, 


.35 
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110816 


12.101 


21 


.673 


25, 


.806 


21.336 


3.515c+14 


0. 


,823 


3839, 


.08 


110817 


12.077 


21 


.437 


25, 


.356 


6.872 


1 1 '■)'■) 1-1/1 

1.133e+14 


0. 


,823 


2631, 


.68 


110824 


10.068 


21 


.454 


25, 


Abb 


A r\ A 

5.404 


5.161e+13 


0. 


,810 


2021, 


A r\ 

.40 


110625 


10.068 


21 


.440 


25, 


.457 


6.877 


6.570e+13 


0. 


,810 


2196, 


.23 


110626 


10.063 


21 


.425 


25, 


.319 


8.846 


O A A r\ 1 "1 o 

8.440e+13 


0. 


,810 


2395, 


.76 


110827 


8.013 


21 


.510 


25, 


.143 


9.017 


4.593e+13 


0. 


,799 


1950, 


.67 


110828 


8.010 


21 


.444 


25, 


.047 


5.885 


2.999e+13 


0. 


,799 


1695, 


.53 


110829 


8.013 


21 


.454 


25, 


.137 


4.908 


2.503e+13 


0. 


,799 


1597, 


.10 


110911 


5.025 


21 


.422 


25, 


.492 


6.604 


1.060e+13 


0. 


,789 


1209, 


.78 


110912 


5.023 


21 


.390 


25, 


.369 


9.779 


1.568C+13 


0. 


,789 


1369 


.93 


110913 


5.023 


21 


.372 


25, 


,387 


12.240 


1.964C+13 


0. 


,789 


1471, 


.19 


110914 


5.020 


21 


.407 


25, 


.206 


8.205 


1.314e+13 


0. 


,789 


1295, 


.76 
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